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Abstract 

We consider a nonlinear oscillator with fractional derivative of the order a. Perturbed by a 
periodic force, the system exhibits chaotic motion called fractional chaotic attractor (FCA). The 
FCA is compared to the "regular" chaotic attractor. The properties of the FCA are discussed and 
the "pseudochaotic" case is demonstrated. 

PACS numbers: 05.40.-a, 05.60.-k, 05.40.Fb 
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I. INTRODUCTION 



It became evident now that the random dynamics can appear as an intrinsic property of 
real systems and that different kind of randomness represents a level of complexity of motion. 
While the Hamiltonian chaos is relevant to the Hamiltonian systems, chaotic properties of the 
dissipative dynamics are very different revealing the chaotic attractors, strange non-chaotic 



attractors, quasi-attractors, etc. 



3[. In this paper we would like to describe one more 
way of the occurrence of chaotic or pseudo-chaotic attractors in dissipative systems. Namely, 
the forced system is a fractional nonlinear oscillator (FNO), i.e. a nonlinear oscillator with 
a fractional derivative of the order < a < 2 with respect to time instead of the second 
order derivative. 

Although the fractional calculus has more than few-hundred-year story, its application to 
the contemporary physics is very recent and, mainly, it is related to the complexity of the 



media in c 



dnetics 



assical and quantum treatment. Let us mention only a few of them: fractional 



6| , wave propagation in a media with fractal properties 



10( 1 , quantum mechanics 



ll| , quantum field theory 



7 



8 



nonlinear optics 



12| . and many others. The forma 



re. 



15 



ated to the contemporary fractional calculus are well reflected in the monographs 



16 



issues 



13 



171. There are different possibilities for interpretation of fractional derivatives. Let 

lall9l|20| 



us mention the probabilistic 



of the considered system 



interpretation 



and the relation to a dissipation 



2lL l22l . |23|, |2J|. The latter one will be related to our paper. 
The paper contains some necessary definitions in Sec. [H] and in two appendices. In 
Sec. mil we demonstrate how the fractional derivative can be related to a dissipation in the 
system. The Sec. llVl is devoted to the main object of the paper: fractional chaotic attractor 
(FCA) and some of its features. We speculate on the existence of fractional "pseudochaotic" 



attractor, i.e. dissipative random dynamics with zero Lyapunov exponent. 



II. DEFINITIONS 

In this section we put some necessary definitions. The left and right tZ)^ Riemann- 
Liouville fractional derivatives of order a are defined as 



1 /•* 



r(n — a) dt^ Ja 



{t - Tr-"-'f{T)dT 



T{n — a) dt 

{n — 1 < a < n) . 



- I {r-tr—^f{T)dT 



(1) 



To construct a solution for a process described by an equation with fractional derivatives, 
one needs the initial conditions that can be written as 



lim„Dr'7W = afe, (fc = l,2,...,n) 



(2) 



or 



liniiDrVW =&fc , (A; = l,2,...,n) 



t-*b 



(3) 



These conditions may have no physical meaning (for a detailed discussion see [16|,|25[). The 
important feature of the derivatives (Q) is that they have no symmetry with respect to the 
time refiection t —t. 



In the following we use the so-called Caputo derivative 



26 



27| defined as 



1 



T{a — n) Ja 



(r-t)"-"-V(")(r)rfr 



r(a — 71) Jt 

{n — 1 < a <n) 



(4) 



with the regular type of initial conditions 



/(^)(a) = au or /^(fo) = 6, ; = 0, 1, {n - 1) 



(5) 



For the left Caputo derivative we use a notation 



(6) 



Left fractional oscillator will be described by the equation: 



QD'^Xi{t) - \xi{t) = 0, 



(7) 



with initial 
transform 



conditions Solution of ((Tj) can be found with the help of the Laplace 



where 



It follows from 



s''X,{s)-\X^{s) = Y.akS 

k=l 



Xi(s) = / e-''xi{t)dt . 
Jo 



k-1 



k=l 



- A 



and the inverse Laplace transform gives 



(8) 



(9) 



(10) 



Mt) = Y.^kr-''E^^^.k+i{\n , (t > 0) 

fc=l 



where 



(11) 



(12) 



is the two-parameter Mittag-Leffler function 
of the system from the present to the future. 

The right fractional oscillator is given by the equation 



28l | . This equation describes a causal evolution 



tD^X2{t)-Xx2{t) = 0, 



(13) 



with the initial conditions Q- The Laplace transform in this case 

X2{p) = r eP'x2{t)dt , (t < 0) (14) 

J — oo 



gives 



n 



p"X,{p)-\X,{p) = Y.auP^-\ (15) 

k=l 

with the corresponding anti-causal solution 

n 

X2{t) = E«A.(-t)"-'^a,a-m[A(-t)"] , (t < 0) . (16) 
k=l 

Similar results can be given for the equations with left and right Caputo derivatives. 
III. DECAY RATE ANALYSIS 



In this section we analyze the dynamics of a nonlinear fractional oscillator (see also 
Let us start from the linear fractional oscillation satisfying the equation 

L)"a; + x = 0, (l<a<2), (17) 

where D" is the Caputo left fractional derivative (see ©)• Equation fll7|) has a solution in 
a form of one-parameter Mittag-Leffler function 

oo / -| \kj.ak 

x{t) = AE^i-n ^ I = AE^,,{-n (18) 

fc=o + 1) 



22 



provided that x'(0) = and x(0) = A (see Appendix 1 for details). According to 
the Mittag-Leffler function may be decomposed into two terms 

Ea.{-n = fa{t)+ga{t), (19) 

where 

^"^ ^"7r7o r2" + 2r°cos(7ra) + l 



gjt) = ^e*™^("/") COS 

a 



t sinf — ) 

The first term is determined by a cut on the complex plane for Mittag-Leffler function, and 
the second one is related to the poles. For a = 2 — £:(£:^l)we can use an expansion over 



fait) ^ -e e 
Jo 



°° r dr 



r2 + l)2 

= --{l-t[ci(t)sint-si(t)cost]}, e<l (21) 
where ci(t) and si(t) are sine and cosine integral functions respectively. For large t the 



function Ei—f^) has algebraic asymptotics of faif) 



22 



28 



f2-e{t) - J. , t>r = — In-, £<1, (22) 
1 — Ij vre e 

i.e. the first term in (fTTHl decays algebraically in time while the second one decays exponen- 
tially. The decay rate of Qait) is defined by its amphtude 

2 

max^„(t) = -e"^*, (23) 
a 

where 

7 = — cos(7r/Q;) . (24) 

It is useful to remark that < 7 < 1 for 1 < a < 2, and 7 ~ ire /A for e <C 1. For a 
comparison the function El—t"') and its envelope ()23|) are presented in Fig. ^ 

The analysis can be extended to a nonlinear fractional oscillator since for small e one can 
apply the averaging over fast oscillations. 

As an example, consider the fractional Duffing equation 

D'^y -y + ay^ = 0, 1 < a < 2, (25) 

where a > is a constant. The steady states are: y = 0, (unstable) and y = ±l/y/a (stable). 




100 



FIG. 1: Rate decay of the linear fractional oscillation with initial conditions x(0) = 1, x'(0) = 
and the index a = 1.95, /i.95(0) = —0.0256. Solid line corresponds to the solution H18() . dash-line 
- to the approximation 

Consider Eq. ()25|1 near a stable fixed point by the change y — >■ w + 1/ ^Ja. Then 



(26) 



Close to the stable fixed point we have a linear equation 



D'^wq + 2wo = 



(27) 



with a solution 



(28) 



and S is a constant. For a = 2 — e and e -C 1 expression ()28|) is well approximated by the 
relation 



Wo{t) 



2B 

2-e 



e-^^*cos(V2t) + 0(e^ 



(29) 



1.4 
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FIG. 2: Rate decay of nonlinear oscillations described by the fractional Duffing equation with 



initial conditions x(0) = 0.5, x'(0) = and the index a = 1.95. The solid line is a numerical 
solution of (|25|) . the dashed line is the envelope curve according to (|29j) . 



with 7 from ()24p . Fig. |21 shows a numerical simulation of Eq. ()25|) in comparison with the 
amplitude obtained from (j^^ . The numerical analysis of Eq. ()25|) is based on the algorithms 



described in (31 1. 

When a = 2, Eqs. (j25j) and (j26|) become undamped. The leading term of the frequency 



of the oscillation w is cuq = 2^/". From 



32j . a nonlinear correction to this frequency is 



3aB^ 



<0, 



(30) 



when aB^ ^ 1, i.e. \uJi\ ^ u. 

From Eq. ()29|). we can present w in the form 



w{t) ^ wi{t) + 6w{t) = Be"^^* cosf(a;o + uji)t + /S] + 5w{t) , 



(31) 
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where (3 is constant and Sw{t) is a correction to u!i{t) due to the term fait) in (j221) that 
describes the polynomial decay of oscillations for fairly large t. 

Concluding this section one state that the fractional generalization of the considered 
nonlinear oscillations is reduced to some effective decay of the oscillations similarly to the 
fractional linear oscillator. The rates of the decay can be estimated and, roughly speaking, 
the larger is the deviation e = 2 — a, the stronger is the decay. 

IV. FRACTIONAL CHAOTIC ATTRACTOR (FCA) 

It is well known that periodic force applied to a nonlinear oscillator, for general situ- 
ation, leads to the Hamiltonian chaotic dynamics in some part of phase space, while the 
same problem with dissipation can lead to the chaotic attractor. One can expect that the 
periodically perturbed fractional nonlinear oscillator should display a kind of chaos that we 
call FCA. Description of this phenomenon is the subject of this section. The basic equation 
is 

D^'x-x + x^ = Fsm{ut), 1 < a < 2, (32) 

where F and u are parameters of the perturbation. For fairly small e = 2 — a, let us 
introduce a cojoint equation 

X + aix — X + x"^ = Fi sin(z/it) . (33) 

Simulation of Eq. (jH^ shows a typical chaotic attractor in Fig. El (a). One can recalculate 
the value of a into a corresponding ai. As a result, in Fig.|ni(b) we see the map of the FCA. 
These figures display a structural difference between CA and FCA. We can assume that 
for small e and up to the terms of e"^ the behavior of the FCA is similar to the CA of the 



10 




X 



FIG. 3: Poincare maps: (a) for the ordinary Duffing equation driven by periodic force: oi = 0.1172, 



Fi = 0.279, vi = 1.0149; (b) for the fractional Duffing equation with the index a = 1.9, = 1, 



F = 0.3, tmax = IStt, driven by periodic force. Number of trajectories in (b) is 2000. 
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FIG. 4: Exponential growth of distance |Aa;| of initially close trajectories {Ax =10 ^, t = 0) for 



the FCA. F = 0.3; i/ = 1. 
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FIG. 5: Reduction of structures for the FCA with the growth of e = 2 — a: (a) a = 1.82, v = 0.985, 
F = 0.39; (b) a = 1.7, u = 0.955, F = 0.6. Number of trajectories is 200. 

cojoint equation. Nevertheless, it seems that this similarity is not complete and some strong 
difference is indicated below. 

The computational time for Eq. (|32j) is very large in order to keep a reasonable accuracy. 
To reduce this time one can consider a short time of computation but with large number of 
trajectories . That is a way how the Fig. EKb) was obtained. It means that the Poincare 
maps for CA and FCA were obtained in different ways. To check the presence of chaotic 
dynamics, we consider growth of the distance |Aa;| for two initially close trajectories for 
Eq. (j32|) . This result is in Fig. 0] and it conforms the presence of a positive Lyapunov 
exponent. Nevertheless, we can see almost regular returns of |Ax| to the initial value 10~^. 

As it was mentioned above, increase of the "fractionality" of the the derivative order 
e = 2 — a can increase the dissipation. This leads to some "reduction of structures" in the 
FCA (see Fig. similarly to the case of CA 2 1 . 



Finally, we present the case of a "dying attractor'' 



2] in Fig. ini Although the phase por- 



trait looks very regularly, the dispersion of initially close trajectory shows randomness with 
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FIG. 6: Dying FCA: a = 1.82, F = 0.3, v = I; (a) phase plane of the FCA; (b) dispersion of 
initially close two trajectories. 

the most probably zero Lyapunov exponent. This case can be related to the fractional pseu- 
dochaotic attractor (FPCA) that has a counterpart (pseudochaos) in Hamiltonian dynamics 
33[. Two comments support this hypothesis: (a) there is no separation of the initially 
close trajectories for fairly large time after which the distance jumped to the order one (see 
Fig. ini(b)); (b) there is no structure of the attractor in Fig. ini(a) even when we strongly 
increase a resolution of the Poincare section plot, i.e. the dimension of the set in Fig. |ni(a) 
is rather one than larger of one as in Fig. El (a). 



V. CONCLUSIONS 

Since the fractional derivatives are time-directed, the equations with fractional derivatives 
slightly different from the integer ones by £ = 2— a <^ 1 can be fairly easy interpreted through 
the regular equations with a dissipation. As a result, fractional nonlinear oscillator behaves 
like the stochastic attractor in phase space, being periodically perturbed. The role of the 
polynomial dissipation is still elusive. It seems that this term leads to a degradation of the 
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FIG. 7: Errors between the numerical calculations with time steps 0.01 and 0.001 (dots); with time 
steps 0.001 and 0.1 (pluses) . 

FCA structure. The error has fast increase with time creating difficulties in the simulations 
(see Fig. [7j). Due to that we can not provide explicit features of the difference between CA 
and FCA. 

The resonant case for the linear oscillator can be interpreted in the way similar to the 
integer derivative case with a = 2 321 (see Appendix 2). Our simulations show that applying 
the fractional calculus, one can gain a compact formulations of dynamics with new properties 
governed by a complexity of the media. 
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Appendix 1 

Here we consider a presentation of the solution to equation (fTTj) . Two independent 
solutions are 

ea{t) = , la{t) = r/'S.,l+„/2(-t") . (34) 

Let u{t) be a function and L.T.u{t) be its Laplace transform, namely 

L.T.u{t)=u{s)= e-''u{t)dt. (35) 
Jo 

Then the Laplace transforms of ea{t) and ia{t) become 

ga-l g-l+a/2 

L.T.eJt) = ; L.T.iJt) = . (36) 

Let us deform the Bromwich path of integration into the equivalent Hankel path. Then 
the loop starts from minus infinity along the lower side of negative real axis, encircles 
|s| = 1 counter-clockwise and ends at minus infinity along the upper side of the negative real 
axis. This allows one to decompose the functions ea{t) and ia{t) into two terms. The first 
contribution arises from two borders of the cut along the negative real semi-axis. The second 
contribution is determined by residues in the poles Sq = exp{jn/a) and Si = exp(— jTr/a). 
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We arrive to the expressions (fT ^ - fpUj) '22"] for the function ea(t). Similarly, we can present 



where 



_ ^-i+a/2(i _ sin(7ra/2) dr 
IT Jo r^" + 2r" cosfvra) + 1 ' 



qa{t) = -e' 
a 



_ ^tcos{n/a) 

For the detailed analysis of these expressions see 



t sm ( — 



(37) 



(38) 



Appendix 2 Driven linear fractional oscillator 



In this section we discuss briefly some features of linear fractional oscillator perturbed by 
an analog of the resonant external force (see also |3^). "Free" and "forced" oscillations of 
the fractional oscillator depend on the index a. The main conclusion is that the dynamical 
response of the driven fractional oscillator is bounded in amplitude for any relation between 
the oscillator frequency and the frequency of the perturbation. The finiteness of the response 
indicates a damped character of fractional derivative. In j2j] the linear fractional oscillator 
is interpreted as an ensemble average over harmonic oscillations because of the interaction 
of the fractional system with the random environment. The intrinsic absorption of the 
fractional oscillator results from the response of each harmonic oscillator being compensated 
by an antiphase response of another harmonic oscillator (see details in 



). This shows the 

main difference between the resonant phenomenon in regular systems and the resonance in 
fractional systems (both linear and nonlinear). 

Let the external force in the linear fractional oscillator equation be described by Mittag- 
LefHer function 



D'^x + u^x = CEa{-ujH") , 



1 < a <2. 



(39) 
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where C is a constant determined by initial conditions. Since the equation is hnear, its 



solution can be presented in the form similar to 



32 



x{t) = AiE^{-ujX) + 



LUn — UJ^ 



2.a\ 



(40) 



where Ai, Ci are new constants. A direct calculation gives 



lim x{t) = AiE^i-u^r) + Ci—E^ ai-uJnt") ■ 



(41) 



For t ^ the second term tends to zero. For l<a<2,t:^lwe have 



2j.a\ 



and the same is for the first term of (^1]) 



^o'r(i - a) 



co^nTil - a) ' 



t > 1. 



(42) 



(43) 



Thus, the secular term is absent for 1 < a < 2. 

The resonance is observed only for a = 2, when Eq. (PT|) is transformed into the usual 
forced harmonic oscillator. Then, the secular term appears 



lim Xa=2{t) = AiE2{-ujie) + CieE2,2{-ojlt^^ 

Cit . 

= Ai COS uot H sm uot . 

2ujQ 



(44) 



Concluding, the linear fractional oscillator forced by a Mittag-Leffler oscillation does not 
lo to any resonance for uj uoq for 1 < a <2. 
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